function dickson_charge_flow_extra

clc;
A1=[0 1 0 1 -1;1 0 1 0 1;1 0 0 1 0;0 1 1 0 0]
NullA1_0=null(A1)
NullA1=NullA1_0(2:end-1,1:end)

A2=[0 1 1 0; 0 0 1 1; 1 0 0 0]
NullA2_0=null(A2)
NullA2=NullA2_0(1:end-1,1:end)

A3=[0 1 1 0;-1 0 -1 -1;-1 1 0 -1];
NullA3_0=null(A3)
NullA3=NullA3_0(1:end-1,1:end)




B=null([NullA1 NullA2 NullA3])

ans1_1=NullA1_0*B(1:2,1)
ans1_2=NullA2_0*B(3,1)
ans1_3=NullA3_0*B(4:5,1)
%sum1=[ans1(1,1);sum([ans1(2:end) ans2 ans3 ans4],2)]

ans2_1=NullA1_0*B(1:2,2)
ans2_2=NullA2_0*B(3,2)
ans2_3=NullA3_0*B(4:5,2)
%sum2=[ans1(1,1);sum([ans1(2:end) ans2 ans3 ans4],2)]



ans1_1+ans2_1
ans1_2+ans2_2
ans1_3+ans2_3
